Nonequilibrium open quantum systems with multiple bosonic and fermionic environments: A hierarchical equations of motion approach

We present a hierarchical equations of motion approach, which allows numerically exact simulation nonequilibrium transport in general open quantum systems involving multiple macroscopic bosonic and fermionic environments. The performance the method is demonstrated for model nanosystem, involves interacting electronic vibrational degrees freedom coupled to baths. results show intricate interplay this scenario both voltage thermally driven processes. Furthermore, use importance criteria improve efficiency discussed.